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Abstract 

Understanding how biological homochirality emerged remains a challenge for the researchers interested in the origin of life. 
During the last decades, stable non-racemic steady states of nonequilibrium chemical systems have been discussed as a possible 
response to this problem. In line with this framework, a description of recycled systems was provided in which stable products 
can be activated back to reactive compounds. The dynamical behaviour of such systems relies on the presence of a source of 
energy, leading to the continuous maintaining of unidirectional reaction loops. A full thermodynamic study of recycled systems, 
composed of microreversible reactions only, is presented here, showing how the energy is transferred and distributed through the 
system, leading to cycle competitions and the stabilization of asymmetric states. 
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Introduction 

The early origin of life cannot be studied without taking 
account the self-organization of chemical networks, 1 the emer- 
gence and antagonism of autocatalytic loops, 2 and the onset of 
energy fluxes driving the whole process. 3 Such chemical net- 
works are especially interesting to understand the appearance 
of biological homochirality — still one of the most prominent 
problems in the origin of life 4 ' 5 — as the destabilization of 
the racemic state resulting from the competition between 
enantiomers and from amplification processes concerning both 
competitors. 6 The onset of biological homochirality can thus 
be conceived as a symmetry breaking phenomenon occurring 
in a nonequilibrium autocatalytic system leading to stable 
non-racemic steady states. 7 

Recent developments in those lines are aiming at the de- 
scription of recycled systems, 8 ' 9 rather than open-flow sys- 
tems. 10 While traditional model systems are totally open, the 
energy input being brought by fluxes of the chiral subunits 
themselves, what we call "recycled systems" are instead closed 
to the chiral subunits but coupled to a flux of energetic 
compounds. In recycled systems, there is a constant number 
of chiral subunits that are driven away from the equilibrium 
state by an active process allowing to transform low-potential 
chiral subunits into high-potential ones. 6 This absence of ad- 
dition/removal process of chiral subunits, which are inherent 
to the Frank's model, implies a necessary interconversion of 
chiral subunits to obtain a change of the initial imbalance 
(see Fig. [T]). Reaction fluxes become thus necessary, going 
through some epimerization reactions, directed from the minor 
configuration towards the major one, that is, in a direction 
opposed to the thermodynamic one. 
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Figure 1: Schematic representation of the deracemization 
process in open-flow and recycled systems. Open-flow systems: 
the autocatalytic construction increases the absolute differ- 
ence between L and D, and maintains the relative difference; 
mutual destruction maintains the absolute difference, and 
increases the relative difference. Recycled system: some of 
minor configuration compounds are converted to the major 
configuration, increasing both the absolute and relative differ- 
ences between L and D. 

While previous studies have focused on the kinetics as- 
pects, 9 the purpose of this work is to understand how such 
chemical dynamic systems can be conceived still obeying 
the thermodynamic laws. Our conceptual framework will 
lay an important emphasis on entropic exchanges. 11 ' 12 This 
approach requires to take into account the microreversibility of 
all reactions, even if they turn out to be quasi-irreversible, and 
to explicit the source of energy. These points can be neglected 
in a kinetic study, 13 but are necessary to compute energetic 
quantities. 

Departing from there, a nonequilibrium Onsager's triangle 
of reactions 14 will be analyzed to point out how chemical 
energy can be transferred into a chemical system and main- 
tain a unidirectional reaction cycle. This will lead to the 



description of nonequilibrium systems in terms of the reaction 
fluxes (that is, the dynamics of the systems), rather than 
the concentrations (that is, their static parameters). Finally, 
the APED system (Activation-Polymerization-Epimerization- 
Depolymerization system 9 ) will be analyzed, showing how a 
recycled system of chiral subunits can be maintained in a non- 
racemic state by consuming chemical energy, and by driving 
this energy towards reaction cycles of autocatalytic inversion 
of chiral centers. An energetic and entropic analysis will show 
how this source of energy can be efficiently consumed. 

General Relations 

Let us consider a system constituted by r chemical reactions 
Rj, involving n compounds Xi. For each compound Xi, n 
fluxes of exchange of matter $i can be established with the 
surrounding: 

$i : ^ Xi (la) 

n k n 

Rr- ,^ »:. v>. (ib) 

i=i 3 i=i 

The exchange flu^ J 15 ^ ipl is positive for incoming compounds, 
negative for outgoing compounds, and zero for non-exchanged 
compounds. 16 The net stoichiometric coefficients are: 

v i,i = v Zj- v i,r ( 2 ) 
If vj- — and vT. 0, Xi is a reactant of the reaction Rj. 
If / and = 0, Xi is a product. If ^~ = vfj =fi 0, 
Xi is a catalyst. If ifj 7^ 0, u^- 7^ and Vi.j 7^ 0, Xi is an 
autocatalyst. If u~- = vf- — 0, Xi does not participate in the 
reaction Rj. 

The rate of appear ance 1 15 ^ Vi in each compound Xi is 
the resultant from exchange fluxes and reaction fluxes (See 
Appendix from Eq. SO] to Eq. [36] for more details): 

r 

= <Pi+J2 Vi 'i'Pj- ( 3 ) 

The purpose of this study is to describe the circulation of 
matter and energy fluxes inside the system. Like appropri- 
ately underlined by Blackmond and Matar, 13 ' 16 all reactions 
must be described as microreversible^^. This means that 
in the equ ilibrium state, each reactions must obey detailed 
balance^" implying that all independent fluxes (pj are zero 
(see the Appendix part "Microreversibility and detailed bal- 
ancing" for more details). As a consequence, direct and 
indirect kinetic constants (kj and fc-j) must be consistently 
related with the thermodynamic equilibrium constant Kj. 

An elegant way to automatically take into account all these 
relationships in complex chemical networks was described by 
Peusner and Mikulecky. 17 An analogy between chemical 
networks and electric networks can be done, in which case 
a chemical reaction can be described by the resistive onset of 
chemical fluxes between reactant and product "nodes". This 
can be done by defining one characteristic parameter per 
element: the Tj for each reaction, and the standard constant 

_ AfG? 

of formation Kf t i — e ht for each compound. In that 
context, a given state of the system is totally specified by 
one variable for each element: the flux tpj of each reaction, 
and the potential Vi of each compound. 

The parameter Vi, characteristic of the state of the com- 
pound Xi is defined by: 

^ = P~ (4) 
= e^r . (5) 



Vi is the activity of compound Xi (taken equal to its concen- 
tration Xi) corrected by its relative stability (reflected by its 
standard constant of formation Kf : i). 

The parameter Tj, characteristic of the reaction Rj is 
defined by (See Appendix from Eq. 27] to Eq. [52] for more 
details): 



(6) 



Fj reflects the general speed of the reaction Rj in both 
directions. The effective constant rate is then equal to Fj 
(high values reflecting fast reaction) divided by the standard 
constants of formation of the reactants (stable reactants 
leading to slow reactions, unstable reactants to fast reactions). 

The expression of the chemical flux for each reaction Rj 
then becomes (Appendix Eq. 1531) : 



<pj = 1 Ik 

i=i 
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(7) 
(8) 
(9) 



This equation is equivalent to a chemical "Ohm's law", de- 
scribing the onset of a chemical flux generated by differences 
of chemical potential between reactants and products. 

The entropy production by one reaction (Appendix Eq. 1581) 

is: 



(10) 
(11) 
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This is equivalent to an energy dissipation, aj being positive 
whatever the actual direction of the flux. The total entropy 
production of the network is: 



The entropy of exchange (Appendix Eq. |62|) is: 

n 

a c = R^^fl InVj. 



(12) 



(13) 



A positive value of a e indicates an increase of internal energy 
Ta c due to the matter exchanges (and symmetrically, a 
negative value indicates a decrease). 

The entropy balance can be computed by: 



Ubai = a — a 



(14) 



The global energy exchange is e — Tabai- When abai > 0, the 
system globally consumes energy. When abal < 0, the system 
globally releases energy. 

A steady state corresponds to: 



fi = -^Vijifi- 

3=1 



(15) 



The exchanges with the surrounding can compensate an 
unbalance of the reactions, maintaining continuous fluxes fj. 
In this case, abal is zero, that is a c — a 1 . The system exactly 



dissipates what it receives, maintaining constant its internal 
state. The maintaining of this state is active, and is different 
from the unique equilibrium state. The matter fluxes f 1 imply 
a continuous creation of entropy during the transformation. 
The matter flux maintains the reaction, allowing it to be 
continuously performed. 

The equilibrium state is a special case of steady state where 
each ipj = and fi = 0, in whic h case the absence of 
exchanges leads to detailed balance^^. There are no exchange 
of any sort with the surrounding, and no entropy is produced: 
<r c , abai and all aj are zero. 

Nonequilibrium Onsager's Triangle 

Coupling a chemical system to a source of energy leads to the 
spontaneous formation of fluxes, possibly leading to positive 
feedback and self-organization. 3 Such couplings are actually 
ubiquitous in all biological systems, where many endergonic 
biochemical reactions are coupled to the hydrolysis of ATP. 
Some nonequilibrium abiotic systems also function on the 
basis of chemical energy transfer. 18 These nonequilibrium 
systems are characterized by a closed system of given chem- 
ical compounds, performing unidirectional loops through the 
consumption of an excess (or clamped) of "fuel" molecule. The 
canonical following example describes how energy can generate 
and maintain an unidirectional cycle of reaction. 

Description of the System. At equilibrium, a reaction loop 
can't be subject to a unidirectional flux of reactions, as all 
reactions become balanced in detail. 14 ' 19 Nonequilibrium 
steady states rely on energy consumption dissipated by en- 
tropy production, leading to cyclic processes. 20 In order to 
understand what happens in these steady state systems, let 
us consider the following system: 
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(16) 



The whole network is represented in fig. [2]A Fluxes of X and 
Y force the reaction from C to A, and thus maintain the 
system {A,B, C} out of equilibrium. The reaction fluxes can 
be expressed as: 



fa = F a (V C V X - VaVy) 

V>1 = Tl (V a - V b ) 

f2 = F 2 (V b - Vc) 

<P3 = F 3 (V c -V a ). 



Steady State. In the steady state, we have: 



(17a) 
(17b) 
(17c) 
(17d) 



Vb = = ipi — (f2 
V x = = tfx - fa 
VY — = f Y + fa 
Va — = f a — fl + f3 



f\=fi = f, (18) 

fX = fa = f B , (19) 

fY = —fa = -f°, (20) 

f3 = f~ f '• (21) 



Note that this comes down to apply the equivalent of the 
Kirchhoff's current law for a steady state chemical system. 21 
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Figure 2: A: Schematic representation of the reaction network 
relative to a nonequilibrium Onsager triangle. B: Steady state 
potential of the system 13a . a, spontaneous transformation 
from A to C; b, activated transformation from C to A; c, 
spontaneous transformation from X to Y, coupled to the 
activated transformation b, the process c forcing the process 
b to compensate for the process a. C: Detailed distribution 
of the entropic fluxes inside the whole network. D: Simplified 
description of the system, emphasizing the coupling between 
the external compounds X,Y and the internal compounds 
A,B,C. 

The relationship between the several fluxes is thus: 

(rr 1 +iY 1 )v? i = V a -Vc (GzEJ + 



(22) 



r 3 -V3 = -(rr 1 + r 2 - 1 )^ i (G5» = 
i i 



(p = 



(23) 

(It2l1)fc(p3|). (24) 



For a positive exchange flux ip c , ip 1 is positive and ip3 is 
negative. There is a net continuous transformation around 
the cycle, initiated from C to A, transmitted from A to B and 
from B to C. ip 1 = <p c + ip3 is the effective flux induced in the 
cycle. if3 is a "leak" of the activation flux. The induced cycle 
flux is equal to the exchange flux (ip' — <p e ) when F3 <C Ti 
and F3 T2, that is, when the kinetics of reaction 3 can be 
neglected compared to the kinetics of reactions 1 and 2. The 
fraction of flux that is induced in the cycle is independent of 
the kinetics of the activation reaction, which only plays a role 
for the intensity of the steady state exchange flux (efficient 
activation reactions imply high exchange fluxes). 

At equilibrium, all the Vi are equal while in the nonequi- 
librium steady state, the exchange of chemical energy allows 
to mainta in V a > V b (Eq. [PTE)) . V b > V c (Eq. IPTc) . V a > V c 
(Eq. ll7dp . and V C V X > V a V y (Eq.[l7aJ) which comes down to: 



and thus V x > V y . Note here that although A is of higher 
potential than C, there is a net conversion from C to A as 
P>a + <fi3 = tp 1 > 0. The transfer of chemical energy allows the 
recycling of low potential compound C back to high potential 
compound A, counteracting the spontaneous evolution in the 
opposite direction (see fig. EJ3). 

Entropic Analysis. The different entropy contributions pro- 
duced or exchanged during the processes are: 
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= R(<p c - p)ln^-. 


(26e) 



The incoming energy, due to matter exchange, is shared by 
the four reactions. These exchanges are represented in fig.[2jl!. 

Simplified Description. The full chemical system can be di- 
vided in two subsystems: the internal system, only subject to 
internal reaction loops, and the exchange system, involving 
compounds that can be exchanged with the surrounding, 
and thus subject to linear chains of reactions connected to 
external sources of matter. These two subsystems can be 
easily determined from the laws of mass conservation. The 
conserved moieties can be computed from the left nullspace of 
the stoichiometric matrix v — [vi,j\ of the reaction network 22 
(the zeros of the stoichiometric matrix are replaced by dots 
for better readability): 
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(27) 



The first vector of the left null space base gives a pool of mass 
conservation for {A, B,C}, and the second vector for {X, Y}. 
X and Y compounds are involved in matter exchanges, so that 
their concentrations are ruled by the external conditions. No 
flux of A, B or C is possible, so the total concentration a + b+c 
is constant (i.e. the system {A, B,C} can be considered as 
closed, in communication with the open system {X, Y}). 

The subsystem {A, B,C} is composed of a single unidirec- 
tional loop of transformations from A to B, then to C, and 
back to A, performed at a constant rate ip 1 (see Fig. [2jD). 
The subsystem {A, B, C} acts as a closed system, maintained 
in a nonequilibrium activity through its coupling with the 
open-flow system {X, Y}. The incoming chemical energy Ta e 
is transferred from {X, Y} to {A,B,C}, and dissipated by 
continuous entropy creation: 



cr a + a\ + <J2 + 03 

e 

(7 . 



(28) 
(29) 



Vc 



(25) 



Provided the system {X, Y} is such that the concentrations 
x and y can be considered as constant (e.g. in the case 
of a large reservoir of energy compared to a small system 
{A, B,C}), the system {A, B, C} is mathematically strictly 
equivalent to a closed system, in which we have first order 
reactions in A, B and C. The corresponding apparent kinetic 
rates would seem not to respect the microreversibility, as they 
implicitly depend on the concentrations x and y; each time 
such theoretical chemical network is built, 23 it is fundamental 
to realize that there is a hidden source of energy in the system. 



Time Evolution. Numerical integration of the several subsets 
of the system were performed using Xppaut, 24 with iff, a = 
1 M, K { ,b = 2 M, Kf, c = 4 M, K f>x = 1 M, K { , Y = 
10 M, ri = 1 M.s" 1 , T 2 = 2 M.s -1 , T 3 = 4 M.s -1 ' and 
T a = 40 M.s" 1 . This corresponds to fci = ki — kz = 1 s _1 , 
k a = 10 M _1 .b _1 , k-i = k-i = 0.5 s _1 , fc_3 = 4 s _1 , 
k- a = 4 M _1 .a -1 , Ki = ^£ = 2, K 2 = = 2 and 

^ = TT q = 10 - Initial conditions were a = 3 M, 6 = 2 M, 
Co = 2 M. x and y = 10 M were maintained constant for 
each case, corresponding to the presence of a large reservoir of 
X and Y connected to the internal system. Then, x was varied 
as a parameter, allowing to tune the quantity of incoming 
energy. 

The evolution of the different fluxes of transformation (pi, 
if>2 and tp' 3 — (p 3 + (p a (i.e. the total fluxes of transformation 
from A to B, from B to C, and from C to A) was examined. 
Under these conditions, three different kind of evolution were 
observed (see Fig. [3): 

Equilibrium (x — 0): There is no energy source. The system 
initially dissipates its excess of internal energy, to evolve 
spontaneously towards the equilibrium state (a = 1 M, 
b — 2M,c = 4M) where the entropy production is 
zero (i.e. the system reaches a maximum of entropy). 
As only spontaneous reactions are possible in absence of 
coupling, the reaction runs from A to B and from B to C 
{(pi > 0, ip2 > 0) and from A to C ((p' 3 < 0). When the 
equilibrium is reached, all the fluxes are zero: all reactions 
obey detailed balance. 

Close to equilibrium (x = 0.01): There is a weak source 
of energy. The evolution is very similar, the system 
evolves towards a steady state that is very close to the 
equilibrium state (a = 1.07 M, 6 = 2.01 M, c = 3.92 M). 
The major difference is observed after some time. The 
entropy production does not go to zero, but diminishes 
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until reaching the exchange entropy, so that the incoming 
energy is totally dissipated. The system thus reaches a 
minimum of entropy production. The fluxes also follow a 
similar evolution as previously ((fii > 0, (p2 > 0, (p' 3 < 0), 
except that they now reach a non zero positive flux <p\ (p 3 
reversing its direction when approaching the steady state 

(pi = tp2 — p'z > 0. 

Far from equilibrium (x = 0.5): There is a strong source 
of energy. This time, the entropy production increases 
before decreasing towards the steady state (a = 2.6 M, 
b — 2.4 M, c = 2.0 M). During the whole time, the 
entropy production remains actually quite close to the 
exchange entropy. A unidirectional cycle of reaction is 
almost instantaneously obtained in the system, (pi > 0, 
ip2 > 0, (p' 3 > 0). The far-from-equilibrium system is 
ruled almost only by the energy flux. 

The distribution of the fluxes for different values of x is given 
in Tab.[T] The exchange entropy grows with the availability of 
the energy source X. It can be seen that the majority of the 
energy is dissipated by the activation reaction, only a small 
fraction of energy being transmitted to the cycle (from 1% 
to 16% in the performed experiments). In accordance with 
Eq. [24] it can be checked that the forced circular flux i always 
represents 14% of the XY exchange flux ip°, whatever the 
value of the incoming flux. 

The APED Model 

Energy Diagram 

The APED model describes a system based on the activation 
of monomers (that can typically be amino acids), the poly- 
merization of activated monomers with unactivated monomers 
or polymers, the epimerization of one end-residue of the 
polymers, and the depolymerization of the polymers. 6 ' 9 All 
these reactions take place at the same time, in a system that 
is closed in terms of monomer residues (they never enter or 
leave the system). The activation reaction is coupled to the 
consumption of chemical energy, maintaining the system in 
a nonequilibrium state. It was shown that such a system 
can lead to a stable non-racemic steady state, with the 
approximation that most of the reactions are irreversible. 9 ' 13 
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Figure 3: Time evolution of entropy productions (A, C, E) 
and chemical fluxes (B, D, F) in nonequilibrium Onsager's 
triangles as defined in Eq. [16] A and B: x — M, C and 
D: x = 0.01 M, E and F: x = 0.5 M. See text for the other 
numerical values. 



Figure 4: Thermodynamic diagram of the different compounds 
engaged in the APED system. 

Let us analyze the behavior of the APED system, limited 
to dimers for the sake of simplicity, appropriately taking 



Table 1: Entropy productions, chemical fluxes and compounds concentration at the steady state in nonequilibrium Onsager 
triangle for different incoming fluxes. 
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into account the microreversibility of all the reactions. 16 In 
order to keep the new system close to the irreversible system, 
it is necessary to correctly choose the parameters so that 
the previously irreversible reactions remain at least quasi- 
irreversible; i.e. the previously neglected reverse reactions 
must be slow compared to the direct reactions. As a con- 
sequence, these reactions must have a very negative A r G°. 
This corresponds to big energy differences between monomers 
and activated monomers, between activated monomers and 
dimers, and between dimers and monomers. The diagram of 
Fig. |4] summarizes the energetic profile of the system. 

A simple source of energy is now explicitly introduced to 
show the origin and repartition of fluxes inside the system. 
In so doing, we have to take into account an activating agent 
X that allows the transformation of monomer into activated 
monomer. The activation reaction will thus become L + X ^ 
V , essentially displaced to the right. Additionally, the spon- 
taneous deactivation of the activated monomer will release a 
low potential molecule Y as a waste. This can be done either 
via a direct deactivation back to monomer L* ^ L + Y , or 
an indirect one through polymerization L* + L ^ LL + Y , 
these two reactions being again essentially displaced to the 
right. A huge difference of energy must exist between X and 
Y to guarantee the quasi-irreversibility of the reactions. The 
transfer of energy is globally brought by the transformation of 
X into Y , as a result of the activation/deactivation process. 

Like for the Onsager's triangle, the system can be de- 
composed into an internal subsystem (corresponding to the 
previous system with explicit energy fluxes) and an exchange 
subsystem: 
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The first vector of the left null space base gives a pool of mass 
conservation for {L, D, L* , D* , 2LL, 2LD, 2DL, 2DD}, and 
the second vector for {L* , D* , X, Y}. X and Y compounds 
are the only ones involved in matter exchanges, so that their 
concentrations, as well as the concentrations of L* and D* 
are ruled by the external conditions. The internal subsystem 
corresponds to the system of monomers, with the conservation 
of their total number: 



mass conservation. All compounds are involved in a linear 
chain of reactions connected to external sources (see the 
Appendix part "recycled and open-flow systems" for more 
details). 

As stated in the original articles, 6 ' 9 this kind of activation 
can correspond to the activation of amino acids into NCA (N- 
carboxyanhydride of a-amino acids). Several activating agents 
X can be used. In this context, carbonyldiimidazole (CDI) is 
a commonly used compound 26 that may have been formed in 
prebiotic conditions. 26 Other prebiotically relevant agents are 
also often cited in the literature, like activation by isocyanate 
and nitrogen oxides, 27 or carbon monoxides in conjunction 
with sulfide compounds. 28 When the activation leads to 
NCA formation, the waste compound Y resulting from the 
spontaneous hydrolysis of NCA is CO2. As required, all these 
X compounds are very reactive, while Y is a very stable com- 
pound. 29 One could additionally focus on the phenomenon 
that leads to the presence and maintenance of X compounds 
on prebiotic Earth (like stable nonequilibrium concentrations 
of nitrogen oxides or carbon monoxide in the atmosphere, 30 
or the production of sulfides by geochemical processes 31 ). As 
long as these external phenomena are present, a chemical 
source of energy would be available, keeping the chemical 
system of amino acids in an active nonequilibrium state. 

Thermodynamic and Kinetic Parameters 

The relations will be written as a function of the parameters: 
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(32a) 
(32b) 
(32c) 



Kact depends on the difference of energy between L* and L 
(or, equivalently, between D* and D). Kdim depends on the 
difference of energy between L and LL (or D and DD). And 
K ep i, in turn, depends on the difference of energy between 
LD and DD (or DL and LL). Thus, the situation is totally 
symmetrical with respect to interchanges between L and D 
above. 

The ensemble of kinetic parameters relative to all the 
chemical reactions involved in the system is given in Tab. El 
This way, the system can be totally characterized, while re- 
maining perfectly compatible with thermodynamic and kinetic 
relations, by the following independent parameters: 

• Thermodynamic: K act , Kdim, 7 

• Kinetic: &a, kn, hp, JtD, fcs, Ur, a, {3 



ctot = l + d+r+d* + 2-ll + 2-ld + 2-dl + 2-dd (31) 



• Internal conditions: I, d, I* , d* , 11, Id, dl and dd 



This property is characteristic of recycled systems. For open- 
flow systems, such as the Frank model, there is no internal 



By choosing the parameters so that all the back reactions 
are negligible comparing to the direct reactions, the system 



Table 2: Kinetic and thermodynamic parameters of the reactions of the APED system. Kinetic parameters of reactions involving 
X or Y compounds are apparent kinetic rates, assuming constant concentration of x and y. The exact symmetrical of each of 
these reactions is also included in the whole network, with exactly the same parameters. 



Reaction 



Activation: 
X + L ^ L 

k-A 

Deactivation: 
L* % L + Y 

k-H 

Polymerization: 

kp.i 

L* +L ^ LL + Y 

fc -P,l 
kp 2 

L* + D ^ LD + Y 

k -P,2 

Depolymerization: 

kD,l 

LL ^ L + L 

k-D,l 
kn,2 

LD ^ L + D 

k-D.2 

Epimerization: 



DL 



LL 



Racemization: 
L ^ D 
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described in the previous study, not taking into account 
microreversible reactions, 9 would be a correct approximation 
of the complete reversible framework described here. We thus 
need, in accordance with the first assumption of energy levels: 



Analysis of the System 

Simulations were performed 24 with c to t = 2 M , kA = 1 s" 



Vx > 



Kact 

K act ■ Vy < K dim < 1. 



(33) 
(34) 
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This corresponds to a big difference of energy between mono- 
mers and dimers, a yet larger difference between monomers 
and activated monomers, and a very high potential of the 
activating agent X. Assuming that x and y are constant 
— maintained by external phenomena, or present in a large 
excess — the reactions involving X and Y compounds corre- 
spond to pseudo first order reactions in monomers or dimers. 
Mathematically, the system is equivalent to a closed system 
in monomer derivatives, maintained in a nonequilibrium state 
by a continuous flux of chemical energy brought by the 
spontaneous reaction from X to Y, these compounds being 
maintained by other external processes. 

The system is subject to a continuous flux of chemicals: 



= R(fx InVx + <py m Vy) • 



(35) 



In the steady state, the incoming flux of X is equal to the 
outgoing flux of Y, so that: 

ipx = tp y = <P e - (36) 
It will be compensated by the consumption rate of X: 

vx = 0, (37) 
ip e = k A {l + d) - k- A (l* +«f). (38) 

As a consequence, the steady-state flux of entropy of exchange 
is: 

a c = Rk A ((l + d)-—±—(l* + d*))ln^. (39) 



KactVx 



v v 



was kept as a variable, allowing to tune the distance from the 
equilibrium. 

Early evolution. Fig. OJ represents the early evolution of the 
system, reaching quickly the racemic steady state, from an 
initial state composed exclusively by a quasi-racemic mixture 
of monomers (initial enantiomeric excess 10~ 5 ). In this short 
period, the system reaches a racemic steady state, i.e. a steady 
state that is issued from the thermodynamic equilibrium 
branch. As observed in the nonequilibrium Onsager's triangle, 
different behaviors are characteristic of equilibrium close-to- 
equilibrium and far-from-equilibrium conditions: 

Fig. [5]A. : there is no energy source, the system reaches 
the equilibrium state. The entropy production decreases 
towards zero. This system dissipates its initial excess of 
energy before reaching the isolated equilibrium state. 

Fig. 03 : there is a weak energy source so that the system 
remains close to the equilibrium state. The evolution is 
very similar to the preceding one, except that the entropy 
production decreases until a minimal value. The system 
continuously dissipates the incoming energy. Similarly 
to the previous system, the system initially dissipates its 
excess of energy. 

Fig. [5]C : there is a strong energy source, the system is far 
from the equilibrium state. The balance of entropy is 
now positive: the system stocks energy — rather than 
dissipating it as previously observed — by acquiring high 
concentration of high potential compounds. 

Following a longer period, we observe a transition of state 
for some values of x (see Fig. [6]): the enantiomeric excess 
abruptly changes to a non-zero value. The system switches 



from the unstable racemic branch towards the stable non- 
racemic branch. This transition is followed by a decrease of 
entropy production . A s imilar dynamical behavior is observed 
in the Frank modeP^' while the chemistry underlying both 
models is rather different. If variations of entropy production 
are typically observed during transiti on b etween state s of 
different stability, as in autocatalyticP^ or oscillating 12 ^ 
systems, the sign of this variation is not necessarily negative. 

Complete Bifurcation Pattern. The same system was studied 
in detail for different values of x, ranging from 10" 3 M to 
10 3 M. The particulars of the different time courses obtained 
are reported in Fig. [Jj A different bifurcation pattern than 
the one usually found is now observed: 

• For low values of x, the system remains racemic. As 
expected, this corresponds to a state where the system 
remains close to the equilibrium since not enough energy 
is available. 

• Then, starting from a critical value of x = 0.0105 M, 
a bifurcation point is observed. The system switches to 
a non racemic branch, a far-from-equilibrium state that 
becomes more stable than the close-to-equilibrium one. 
This corresponds to the expected behavior happening 
when enough energy is available. 
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Figure 6: Decrease of the entropy production in an APED 
system during the transition from an unstable racemic steady 
state towards a stable non-racemic steady state, in the case 
x = 0.5 M. 

• However, for the value x = 0.15 M, the enantiomeric 
excess actually reaches a maximum, and then decreases 
back to a second bifurcation point at x = 2.11 M. The 
system then switches back to the racemic branch as the 
only remaining stable steady state. This behavior can 
actually be also observed in the com plete Frank model 32 
or in open-flow autocatalytic system a 12 a l 

Variation of Entropy Production during the Transitions. A 

monotonous increase of the incoming energy Ta° as a function 
of x is observed. In the non racemic stable state — when it 
exists — the system consumes up to 2% less energy than on 
the racemic branch. This corresponds to the decrease of the 
total entropy production during the transition between steady- 
states. This decrease is mostly due to a decrease of entropy 
production by activation/deactivation reactions, and by the 
epimerization reactions. These decreases are small in absolute 
value, but the relative variation of entropy production is very 
important for the epimerization reactions (up to 60%). 

We thus observe a small change on the global entropy 
production of the system, that allows a large local change 
centered in the epimerization reactions — the engine of the 
deracemization process, as they are the only reactions able to 
change the L/D balance in a system in which the total number 
of monomers is constant — and a local increase of the entropy 
production by the polymerization/depolymerization processes 
— that is, a redirection a larger part of the energy fluxes 
towards the autocatalytic loops. Globally, the switch to the 
far-from-equilibrium branch allows the use of the incoming 
energy to both produce more dimers and to spontaneously 
favour the stable non-racemic state. 

Evolution of the Entropy Production with the Energy Flux. 

The global entropy production follows the global energy con- 
sumption and thus also monotonously increases as a function 
of x. For all systems, the energy is mostly dissipated by 
the activation and deactivation reactions. More interestingly, 
the entropies produced by the polymerization, depolymer- 
ization and epimerization reactions are actually increasing 
when stable non-racemic steady states exist, but decrease 
to zero after the second bifurcation point. That is, after 
the second bifurcation point, the activation/polymerization/ 
epimerization/depolymerization loops disappear, giving place 
to only activation/deactivation loops. 

Interpretation. When more and more energy is introduced 
into the system, more and more compounds are maintained 
at high level of energy. As the key compounds (i.e.: the 
dimers, since they allow differences of energy between the 
diastereoisomers) are at intermediate energy levels, the der- 
acemization will only occur for an interval of energy flux. Not 




Figure 7: Enantiomeric excess (ee), diastereoisomeric excess of dimers (eed — ^ ^[ll] + [dd] + \ld] + [dl] ^ ) an< ^ entropy production in 
the steady states as a function of the incoming flux of energy in APED systems (see text for parameters). Solid curves correspond 
to values for stable steady states, and dotted curves to unstable steady states. A, H, P, D, E and R stand for activation, 
deactivation, polymerization, depolymerization, epimerization and racemization, respectively. The gray area corresponds to the 
range where non-racemic states are stable. 



enough energy can not produce enough dimers while too much 
energy destroys the dimers in favor of activated monomers (see 
Pig.®. 

Before the first bifurcation, the system is very close to 
equilibrium, the most populated states being the lower energy 
states (i.e. the free monomers). There is just not enough 
energy to be directed towards the autocatalytic cycles. 

After the second bifurcation point, the situation is reversed. 
The high flux of energy is almost only directed to the activa- 
tion/deactivation loops. The most populated states are the 
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Figure 8: Steady-state composition in APED systems for 
different energy fluxes. Top figure: steady state concentrations 
of free, polymerized and activated monomers as a function of 
the concentration of the activating agent X. Bottom figures: 
schematic representation of major cycle fluxes before the first 
bifurcation point, between the two bifurcation points, and 
after the second bifurcation point; the letter sizes represent 
the relative quantity of compounds; the arrows represent the 
major cycle fluxes. The gray area corresponds to the range 
where non-racemic states are stable. 



highest energy states. This is characterized by a saturation 
phenomenon with most of the monomers being maintained in 
the activated form, rather than in dimers or free monomers. 
Just a small amount of energy is directed towards the auto- 
catalytic loops. This phenomenon can be compared to what 
is observed in open-flow Frank-like systems. 32 When matter 
fluxes are too important, the system remains racemic: the high 
energy incoming compounds have no time to react before being 
flushed outside the system. Similarly to what is observed here, 
the excess of energy flux (due to matter exchanges) maintains 
high concentration of high energy compounds, preventing 
them to react in autocatalytic reactions. 

It is only between the two bifurcation points that the 
autocatalytic loops can be efficiently performed. A substantial 
quantity of dimers is maintained, so that the epimerization 
reaction can effectively take place. The energy is thus 
efficiently used to generate a non-racemic steady state. 

Conclusion 

The full thermodynamic study of the APED system empha- 
sizes the great importance of the energy transfers. Energy 
must flow into the system in order to stabilize the non racemic 
state. But just the consumption of energy is not sufficient: 
energy must be efficiently used. It must be directed towards 
the interesting elements of the reaction network, that is the 
autocatalytic loops, able to induce positive feedbacks inside 
the system. This point has to be focused on when designing 
chemical experimental systems that can present such non- 
linear behaviors. 

This redirection of some part of the free energy available 
in the environment (here, the difference of chemical potentials 
maintained between X and Y compounds) towards some inter- 
nal "mechanisms" (here, the chiral autocatalytic cycles) is one 
of the key conditions allowing the emergence of nonequilibrium 
properties (here, a non-racemic steady state). Studies of the 
origin o f life h ave to take into account such type of "energy 
focus" 33 | 3 a | 3 b | . Real chemical systems are generally away 



from the equilibrium state as they continuously communicate 
with their surrounding. Several abiotic systems can be found 
in very different environments (interstellar space, 34 planet 
crust 31 or atmosphere 30 ) and diverse forms of available free 
energy can be found in these environments: by means of 
these external energy sources, nonequilibrium concentrations 
of activated compounds can be maintained. 

In this context, the challenge is to better understand how 
this available energy can be passed down to other chemical 
systems and used so that it leads to self-organization, in a sim- 
ilar way as it is performed by metabolisms. This corresponds 
to understanding the emergence of protometabolic systems. 35 
When a chemical system is connected to chemical sources of 
energy, internal reaction loops are possibly maintained, poten- 
tially leading to the emergence of autocatalytic behaviors^. 
Such nonequilibrium chemical systems are prone to self- 
organization and self-maintenance. 1 ' 2 Thus, the transition of 
such protometabolic systems towards replicative systems can 
be seen as a process of spontaneous self-organization under 
energetic pressures. Rather than opposing self-organization 
and natural selection, 36 it may help to understand how 
natural selection emerges from self-organization. 37 Avoiding 
the vitalist idea of biochemical systems behaving in a funda- 
mentally different way than abiotic chemical systems, it leads 
to the description of general chemical systems under similar 
principles. 

As a consequence, if a given system permit the contin- 
uous activation of amino acids with a sufficient energetic 
efficiency, 26 polymerization/depolymerization cycles will nat- 
urally be established. This would correspond to a very simple 
self-organized protometabolism, able to consume energy, to 
generate some complexity by synthesizing and maintaining a 
dynamic pool of oligopeptides, and potentially being able to 
evolve towards a stable non racemic state. If such scenario 
could have been realized, it would imply mirror symme- 
try breaking to be one natural property of self-organizing 
protometabolisms. Of course, strong precautions must be 
taken when dealing with such complex chemical systems, 
as emphasized by Orgel. 38 The robustness, stability, and 
chemical relevance of this scenario must be thoroughly tested 
and examined. 

Acknowledgment: K. Asakura, R. Pascal, D.K. Kondepudi and 
A. Brandenburg are gratefully acknowledged for comments and 
discussions. 

Appendix 

Microreversibility and Detailed Balancing. It is important to 
understand that chemical reactions are reversible at the microscopic 
level, and irreversible at the macroscopic level. The microscopic 
reversibility (or microreversibility) is a consequence of the time- 
reversal symmetry of the fundamental laws of physics: for a 
given molecule subjected to a chemical transformation, the inverse 
transformation can be obtained by following the exact pathway but 
in the opposite direction. The detailed balancing is a macroscopic 
property that can be deduced from microreversibility: when the 
system is at thermodynamic equilibrium, reaction rates in both 
directions cancel for each reaction (implying the relationship given 
in Eq. [6j. This transition from a microscopic to a macroscopic 
description — involving statistic treatment and the loss of infor- 
mation about individual molecule behaviors — explains how as 
chemical system becomes irreversible at the macroscopic level, in 
the sense that chemical reactions always evolve spontaneously from 
a nonequilibrium state towards an equilibrium state, never in the 
opposite direction. More details on this subject can be found in the 
article of Boyd. 19 

In that context, there is thus no detailed balancing in a nonequi- 
librium system, in the sense that not all the individual fluxes have 
to be cancelled: there can be the onset of chemical fluxes inside the 
system. 20 However, this does not imply that there is a breakdown 
of time-reversal symmetry at the microscopic level (as it may be the 
case in some specific systems 39 ). In a nonequilibrium system, there 



can be chemical fluxes, but kinetic constants remain linked through 
the thermodynamic equilibrium constants. 

Kinetic Relationships. The system exchanges matter with its sur- 
rounding. The variation in concentration due to external exchanges, 
for each incoming or outgoing compound, is: 



d c X{ 

dt 



(40) 



Note that <^>| can be positive or negative, depending on whether 
the component comes into the system or leaves it. Internal trans- 
formations occur through the chemical processes. The variation in 
concentration due to the internal transformations is: 
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The total variation in the concentration of each chemical species is 
thus given by: 

(45) 



dt 

dt dt 



(46) 



Microreversibility. At equilibrium, the microreversibility implies 
the detailed balancing of the reaction (i.e. the fluxes are canceling 
in both directions of each reaction): 

<pj = <pj, (47) 
^ = (48) 

K -3 i=l 

= K 3 . (49) 

The equilibrium constant can be expressed as a function of the 
standard constants of formation of the compounds involved in the 
reaction: 



(50) 



It is then possible to express a relation between the parameters of 
the direct reaction with the parameters of the indirect reaction: 



(51) 

= rj. (52) 

This leads to the expression of Fj , that is a characteristic parameter, 
combining kinetic and thermodynamic properties of the microre- 
versible reaction. 

Relation fluxes-potential. 
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The parameters relative to the reaction (S^) are separated from the 

V ■ ■ 

parameters relative to reactants (JJ V i ' 3 ) and from the parameters 

„+. 

relative to the products iJ\V- *' J ). This equation is analogous to 
an "Ohm's law", in which the reaction flux (<Pj), as an analogous to 
an electrical intensity, becomes a function of the potential between 
reactants and products (V; parameters). Similar laws are given 
in the literature, but the relationship is sometimes given as a 
function of chemical potentials fii rather than the Vi parameters, 



in which case this law becomes true only close to equilibrium. 40 
The relation of Eq. [53] is true even in nonequilibrium states, and 
is rather similar to the description of Peusner et al., 17 except that 
the Vi potentials are calculated relatively to the standard state — 
thus the introduction of the Kf i — rather than relatively to one 
arbitrary node. The advantage of this description is to be totally 
symmetric for all compounds, not depending on any numbering, 
and is thus more general. 

Entropic Analysis. The system will continuously dissipate energy, 
through the chemical reactions. The entropy production by unit of 
time (Tj ;, per each reaction taking place in the system, is equal to: 41 
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It can be noted that whatever the sign of <pj, in accordance with 
the second principle, crj is always positive (zero in the limit case 

Vi =o)- 

The total production of entropy by all internal transformations 
is the resultant of all these chemical processes, that is: 



(60) 



At the same time, the system exchanges entropy through matter 
exchange with the surrounding: 
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Recycled and Open-Flow Systems. An open-flow system is 
characterized by the absence of internally conserved compounds. 
All of them are linked to matter exchanges, so that there cannot be 
any mass conservation rule. The following simple Frank model is 
taken as an example: 
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The conserved moieties of this system are: 
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All the compounds are engaged in a linear succession of reactions 
connected to exchange fluxes. 

In order to recycle this system, the products must be activated 
back to the reactants, that is to replace the fluxes of Eg. I63al and 
Eg. l63el bv: 
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Which leads to: 
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This now turns out to be an internal subsystem with mass 
conservation of a + I + d + 2p, connected to the flux of matter 
{X, Y}. However, this system is still not a recycled system of chiral 
subunits, as there is no mass conservation of £ + d, except in the 
case where the A and P compounds can be neglected compared to L 
and D. In this latter case, the process of destruction/construction 
of chiral subunits becomes equivalent to an activated racemization 
process. 
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